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Abstract 

In this short note we discuss ordinary differential equations (ODEs) which lin- 
earize upon one (or more) differentiations. Although the subject is fairly elementary, 
equations of this type arise naturally in the context of integrable systems. 

2000 MSC: 34A05. 

Keywords: Linearization of ODEs. 

1 Introduction 

Let us consider a linear n-th order ODE with the general solution 

u(x) = a\f\{x) + ... + a n f n (x), (1) 

which is a linear superposition of n linearly independent solutions fi(x). Imposing a 
non-linear relation among the coefficients, 

F(m, ...,a n ) = 0, 

one obtains an (n— l)-parameter family of functions u(x) which automatically solve a non- 
linear ODE of the order n — 1. By construction, this ODE linearizes on differentiation. 
Imposing two relations among the coefficients, one obtains an ODE of the order n — 2 which 
linearizes on two differentiations, etc. Using this simple recipe one can generate infinitely 
many examples of linearizable equations. This note was motivated by the observation that 
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equations of this type arise naturally in the context of integrable systems. The paper is 
organised as follows. 

Sect. 2 contains a list of examples of ODEs which linearize on differentiation. These equa- 
tions appear in the construction of exact solutions of integrable PDEs, in the classification 
of integrable hydrodynamic chains, etc. 

In Sect. 3 we derive necessary and sufficient conditions for an ODE to linearize upon a 
finite number of differentiations. 

In Sect. 4 the general form of non-linear ODEs linearizable by a differentiation is discussed. 
It is obtained by imposing non-linear constraints among first integrals of a linear equation. 
Parallels with the theory of linear invariant subspaces of nonlinear differential operators 
are briefly discussed. 

We point out that the problem of linearization of nonlinear ODEs has attracted a lot 
of attention in the literature. The conditions of point and contact linearizability of second 
and third order ODEs were first studied by Lie [13], see also [1] (p. 38), [2] (p. 202), 
[3, 5, 11] and references therein. A short remark on 'integration via differentiation' can be 
found in Kamke [12], sect. 4.14. We emphasize that in this paper we are concerned with a 
different concept of 'linearizability on differentiation'. Notice that while the linearizability 
by point or contact transformations is closely connected with symmetry properties of the 
equation, this is not true in our case: the symmetry group of the original nonlinear equation 
can be trivial, becoming non-trivial for a linear equation obtained on differentiation. 

2 Examples 

Example 1. As shown in [10], the construction of 'follyton' solutions of a nonlinear 
system associated with a fourth order self-adjoint spectral problem, reduces to an ODE 

u "" u _ „'V + I n "2 = I C V 

2 2 

c = const. On differentiation this equation becomes linear, u'"" = c 4 u' , with the general 
solution 

u = ao + a\ sinh cx + 02 cosh cx + 03 sin cx + 04 cos ex. 

The substitution of this ansatz into the equation leads to a single quadratic relation among 
the coefficients, — a\ + a 2 . + a 2 + a\ = \oq. 

Example 2. The classification of integrable Hamiltonian hydrodynamic chains associated 
with the Kupershmidt-Manin bracket reduces, in a particular case, to a solution of the 
nonlinear ODE [6] 

(4c 2 xV - Y2c 2 xu -l-ax + 2cx 2 )u"' + (a - 4cx + Y2c 2 u + 4c 2 xu' - 2c 2 x 2 u")u"+ 

(4c - 8cVK -5 = 0, 

here c, a are arbitrary constants. Remarkably, this complicated equation linearizes on 
differentiation, taking the form 

(4c 2 xV - Ylc 2 xu - 1 - ax + 2cx 2 )u"" = 0. 
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Leaving aside the possibility that the coefficient at u"" equals zero (see [6] for a complete 
analysis), we conclude that u must be a cubic polynomial, 

u = ao + a\x + ci2X 2 + CI4X 3 , 

where the constants satisfy a single relation 12ci4 — 8cai + 16c 2 (a 2 — 3a2ao) — 4a2a + 1 = 0. 

Example 3. Another subclass of integrable hydrodynamic chains from [6] is governed by 
the ODE 

8x 2 u"'u + 8xu"u - Ax 2 u" 2 - u' 2 - I2u = 
which linearizes on differentiation, 



8x 2 u"" + 2Axu" + Gu" - 12 = 0. 
The general solution is given by the formula 

u = x 2 + ao + a±x + a2X 1 ^ 2 + 03a; 3 / 2 

where the constants satisfy a single quadratic relation 12ao + a 2 — 30,20,3 = 0. 
Example 4. One of the versions of equations of associativity [4] reads as 

Looking for solutions in the form F = ^ 3 u(x), x = 77/^, one arrives at the ODE 

Quu" - Axu'u" + 2xu" 2 - 2u'u" = 1, 

which takes the form (Gu — 4xu')u"" = after a differentiation. The case u"" = leads to 
the general solution 

u = ao + a\x + U2X 2 + CL3X 3 

where the constants satisfy a single quadratic relation 9aoa3 — a\02 = 1. In terms of F 
these solutions correspond to polynomials cubic in £ and 77. The case Gu — 4xu' = leads 
to u = cx z l 2 . The corresponding F is given by the formula F = c(^rj) 3 ^ 2 where c = z^p- 

Example 5. The third order ODE, 



„, _ „ (a + l)u' - 2xu" 
~ su ( s + l)((s + 2)u-2xu'y 



s = const, arises in the classification of integrable Hamiltonian hydrodynamic chains asso- 
ciated with Kupershmidt's brackets [7]. It possesses a remarkable property: for parameter 
values s = 1,2,3, ... this equation linearizes on exactly s differentiations. Thus, for s = 1 
the differentiation of (2) implies u"" = 0, so that the general solution is 

u = ao + 3aix + 3a2X 2 + 03a; 3 , 

where the constants satisfy a single quadratic constraint 0,003 — a\02 = 0. For s = 2, 
differentiating (2) twice, we arrive at = with the general solution 

u = oq + Aa\x + Ga2X 2 + 4a3X 3 + 04X 4 , 
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where the constants a, L satisfy a system of quadratic constraints 

ao£i3 — a\a2 = 0, 0,10,4 — 0203 = 0, 0004 — a 2 , = 0; 

notice that these constraints specify a determinantal variety characterized by the require- 
ment that the rank of the matrix 

/ a ai a 2 \ 
y a 2 a 3 a 4 y 

equals one. The mystery of this example is unveiled by the formula for its general solution, 

u = a(x + c) s+2 + b(x - c) s+2 , 

which is valid for any s; here a, b, c are arbitrary constants (we thank A. P. Veselov for this 
observation) . 

3 Necessary and sufficient conditions for the linearizability 

In this section we demonstrate how to derive necessary and sufficient conditions for a 
non-linear ODE to linearize on one (or more) differentiations. The procedure is fairly 
straightforward and can be readily adapted to particular situations. 

3.1 First order ODEs which linearize on one differentiation 

Let us characterize first order equations 

v! = f(x,u) 

which imply a linear equation, 

u" = a(x)u + b(x)u + c(x), 

on one differentiation. Thus, we have f x + f u f = af + bu + c. Differentiating this 
relation twice with respect to u, and introducing F = f x + f u f, one obtains F uu = 
afuu- Differentiating this by u once again one has F uuu = af uuu . Thus, the required 
linearizability condition takes the form 

■Fuuu fuu = fuuuFuu] 

see Sect. 4 for the general (implicit) form of all such right hand sides f(x,u) (formula 
(6))- 

3.2 Second order ODEs which linearize on one differentiation 

Let us characterize second order equations 

u" = f(x,u,p), p = u', 
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which imply a linear equation, 

u" = a(x)u" + b(x)u + c(x)u + k(x), 

on one differentiation. Thus, we have f x + f u p + f p f = af + bp + cu + k. Applying to this 
relation the operators d 2 , d u d p , d 2 , and introducing F = f x + f u p + fpf, one obtains 

■Fuu = ®fuu, Fup = Q>fup, Fpp = afpp, 

or, equivalently, d 2 F = ad 2 f (here the second symmetric differential d 2 is calculated 
with respect to u and p only). Differentiating this once again by u and p one obtains 
d 3 F = ad 3 f. Thus, the required linearizability condition takes the form 

d 3 F d 2 f = d 3 f d 2 F. 

3.3 First order ODEs which linearize on two differentiations 

Here we characterize first order equations 

u = f(x,u) 

which imply a linear equation 

u" = a(x)u" + b(x)u + c(x)u + k(x) 

after two differentiations. Introducing F = f x + f u f and G = F x + F u f, we have G = 
aF+bf +cu+k. Differentiating this twice with respect to u one obtains G uu = aF uu +bf uu . 
This implies G uuu = aF uuu + bf uuu and G uuuu = aF uuuu + bf uuuu . Thus, the required 
condition is 





Guu 


± uu 


fuu 


det 


G U uu 


VP 

1 uuu 
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Guuuu 


VP 

1 uuuu 


fuuuu 



= 0. 



In all of the above examples, the linearizability is characterized by differential relations 
which must be satisfied by the right hand side of the equation. As we demonstrate in 
the next section, these differential equations can be integrated in closed form, leading to 
(implicit) representation for all linearizable equations. 

4 General form of linearizable equations 

All equations linearizable by a differentiation can be obtained by imposing functional rela- 
tions among first integrals of linear equations. Since the first integrals can be parametrized 
explicitly by arbitrary functions of the independent variable x, this provides a general for- 
mula for equations which linearize on differentiation. 

Example 6. Let us describe all ODEs which reduce to u"" = after one or two differen- 
tiations. The basis of first integrals consists of 

I L = u '\ I 2 = xu" - u", I z = x 2 u" - 2xu" + 2u, h = x 3 u" - 3x 2 u" + 6xu - 6u. (3) 
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Any third order equation which reduces to u"" = after one differentiation can be repre- 
sented by a single relation among the first integrals, 

F(h,I 2 ,I 3 ,I 4 )=0. 

Any second order equation which reduces to u"" = after two differentiations can be 
represented in implicit form by two relations, 

F(h,I 2 , 1 3 , h) = 0, G{h, I 2 , h, h) = 0; 

one has to eliminate v!" to obtain the required second order equation. 
In general, let us consider a linear ODE 

L[u] = n (n) + h(x)u( n -V + . . . + b n ^{x)u' + b n (x)u = b(x). (4) 

Let fo(x) be its particular solution, and let fi(x),...,f n (x) be a fundamental system 
of solutions (FSS) of the corresponding homogeneous equation. A complete set of first 
integrals for the equation (4) can be taken in the form 

T r i W[fi, ... , fi-i,U - f , fi+l, ■■■ , fn] ■ -, ,^ 

m = w [fl ,...j n ] ■ , = 1 '-' n ' (5) 

where W[- \ denotes the Wronskian of the functions indicated in square brackets. Indeed, 
for an arbitrary solution f(x) = fo(x) + a\fi(x) + ... + a n f n (x) we have Ii[f(x)] = a,i for 
i = l, ... ,n, i.e., all {/j} take constant values. Applying (5) to the equation u"" = with 
fo(x) = 0, fi(x) = x 1-1 , i = 1, ...,4, one obtains first integrals which coincide with (3) up 
to constant factors. 

For example, in the case n = 2 the general form of first order equations linearizable by 
one differentiation is represented via arbitrary functions fo(x), fi(x), f 2 {x) and F as 

F(Ji[«],/ 2 [«]) = (6) 

where 

T M _ W[u-f J 2 ] _ (u-h)f' 9 -{u-f )'f 2 r n _ W[fi,u-f ] _ -(«-/o)/f + («-/o)7i 

J U U J - W[h,h] ~ hf' 2 -}[h ' J 2L«J - w[h,h] ~ hf' 2 -f[f2 ' 

The differentiation of (6) yields 

~m}th]? 2 ( (u " fo) " w[fi > h] ~ {u ~ fo)'(w[fij2}y + (u - h)w[f[, r 2 i) = o, 

leading to the linear equation 

W[f 1 J 2 ,u-f }=0, or W[f 1 ,f 2 ,u]=W[f 1 J 2 ,f \. 

This construction generalizes to the case of arbitrary n in a straightforward way. In 
particular, equation (6) provides the general (implicit) form for equations discussed in the 
Subsection 3.1. Similar representations can be obtained for all other cases from Section 3. 
The linear span of the functions fi(x), 

W n = L{f 1 (x),...,f n (x)}, (7) 



6 



represents the linear space of solutions to the homogeneous equation L[u] =0 correspond- 
ing to (4). The space W n is said to be invariant with respect to a differential operator F, 
if i^fWn] C W n . A systematic study of operators preserving a given subspace was initiated 
in [8] in the context of constructing explicit solutions for nonlinear evolution equations. 
The general form of operators preserving the subspace (7) is given by 

n 

F[u}=J2A i (h,...,I n )fi(x), (8) 
1=1 

where A l (h, . . . , I n ) are arbitrary functions of the first integrals of the equation L[u] = 
(see [14] and [9] for more details). Given an operator F of the form (8), we introduce the 
equation 

F[u] = 0, (9) 
and look for its solutions in the form (1). The substitution of (1) into (9) yields the identity 

n 

Y,A i (h,...,I n )f i (x)=0 
i=i 

implying that 

A i (I 1 ,...,I n ) = for i = l,...,n. (10) 

This system imposes relations on the coefficients Oj. 

Most of the previous examples fit into this scheme. For instance, the equations from 
examples 2 and 4 are written via the first integrals (3) as 

F[u] = 2C 2 {hh - if) + 2Ch -h-aI 2 -^ = 

and 

F[u] = -I A h + hh -1 = 0, 

respectively, where I are given by (3). In both cases the system (10) reduces to a single 
relation. 

Similarly, for the equation (2) rewritten as 

F s [u] = {s + l)[(s + 2)u"'u - su"u'\ - 2x[(s + l)u'"v! - su" 2 ] = 0, (11) 
one obtains, for s = 1, 

Fi[u] = \{hh-hh) = 0\ 
this again leads to a single relation (10). In the case s = 2 we have a representation 

F 2 [u] = ^[(hh ~ hh)x 2 + (Iil 5 - l!)x + (hh ~ hh)] = (12) 

via the first integrals of the equation = 0. In this case the system (10) 

formally consists of three relations, 

hh ~ hh = 0, hh - II = 0, hh - hh = 0, 
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however, only two of them are functionally independent. 

Remark. In accordance with [15], every operator (8) of the order n — 1 — k, admitted by 
the equation 

u (n) = 0, (13) 

is expressed in terms of the differences 

4 = xjt 1 - Jf+Y, i = l,...,n-k, (14) 

with jf = Ii, i = 1, . . . ,n ( ij are first integrals for (13)). All these expressions are of 
the order n — 1 — k, and satisfy the identity D k+1 jf = on solutions of the equation 
(13). Setting n = s + 3, k = s — 1, we obtain that any third-order operator admitted by 
the equation u( s + 3 ) = is defined via the functions J/ -1 , i = 1, . . . ,4. For instance, the 
operator F s from (11) is represented as 

Fs[u] = \{Jt 1 4- l -jr 1 Ji- 1 )- 
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